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Abstract 

We theoretically study the current-perpendicular-to-plane magnetoresistance of a domain wall 
confined in a current-confined-path (CCP) structure made of a nano-oxide-layer (NOL). In order 
to calculate the MR ratio of the system, the continuity equations for charge and spin currents are 
numerically solved with the three-dimensional CCP geometry by use of finite element method. It 
is confirmed that the MR ratio is enhanced by the CCP structure, which is consistent with the 
experimental results. 
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Current-perpendicular-to-plane giant magnetoresistance (CPP-GMR) has attracted much 
attention for its potential application as a read sensor for high-density magnetic recording. 
In order to realize a high-density magnetic recording, we need MR devices with high MR 
ratio and low resistance area product (RA). Although the RA value of a CPP-GMR system 
is much smaller than that of a tunneling magnetoresistance (TMR) system, the MR ratio 
of a conventional CPP-GMR system still remains a small value of a few %. Much effort has 
been devoted to increasing the MR ratio of the CPP-GMR system. One of the candidates for 
such a device is the CPP-GMR spin- valve with current-confined-paths made of a nano-oxide- 



layer (NOL) with a lot of small metallic channels 



The MR ratio of 10.2% is obtained 



by the GPP spin-valve with non-magnetic Cu nanocontacts. This MR enhancement by the 



.a. 



CCP structure was theoretically studied in detail in Refs. 

Recently, Fuke et al. showed that such MR enhancement appears even in the system with 
ferromagnetic nanocontacts. They obtained the MR ratio of 7-10% by the CoFe nanocon- 

n 

tacts in NOL with CCP structure M- It seems that this MR effect originates from the spin 
accumulation around the domain wall confined in the NOL. 

In order to analyze this experimental result, we calculate the MR ratio of the CCP 
system with a domain wall. The MR ratio is determined from the distribution of the 
spin accumulation around the domain wall, which is obtained by solving the spin diffusion 



equation 



7|] with the use of finite element method 



SB 



Let us begin with a brief introduction to the theory of the spin diffusion equation 
which we solve to determine the MR ratio of the system. According to the Ref. the 
electric current je and the spin current jm can be written as 

je = 2C^E - 2D • Vm, (1) 
L = 2CE - 2/^0 Vm, (2) 

where we choose the unit such as the charge of an electron and the Bohr magneton are set 
tobel: e = /iB = l. Here we omit the term proportional to the derivative of the charge 
accumulation. The vectors denoted by the arrow v represent those in real space and the 
bold font V represent those in the spin space. E is the electric field and m is the spin 
accumulation. The conductivity C and the diffusion constant Z) are written in the spinor 
form as C = Cq + a ■ C, D = Dq + a - Y), where the a are the Pauli matrices. C and D are 
proportional to the spin polarization parameters [3 and /5' as C = /SCqM and D = /S'DqM, 



where M is the unit vector in the direction of the local magnetization. 

In order to derive the equation to determine the spin accumulation m, we eliminate the 
electric field E from the equations ([T]) and ([2]) and obtain 



Vm-/3/5'M(M- Vm 



For convenience we introduce an operator A acting on the spin space as 



(3) 



A = -2Do 



ll \ 



\ V 



-(3(3' 
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(4) 



where we write the components of the local magnetization as M = (Mi, M2, M3). With the 
use of this matrix A, the expression of the spin current is simplified into 



j„ = /?Mje + AVm. 

The continuity equation for the spin accumulation takes the form Q 

dm 



(5) 



m 



dt 



-V-j„^-(J//i)mxM , 



(6) 



where J is the coupling constant of the interaction between the spin accumulation and the 
local moment, and Tgf is the spin-fiip relaxation time of the conduction electron. In order to 
obtain the solution for the stationary state, we set the left hand side of Eq. ([6]) to zero. 



= V • (J/n)m X M + 



m 



m 



V • /3Mje + AVm + {J/h)m x M H . 



(7) 



Introducing the matrix B on the spin space as 

ll \ I 



(J/n) 



M3 -M2 
-M3 Ml 
M2 -Ml 



V V V 

we obtain the equation for the spin accumulation m to be solved 

V • f/?M Je + AVm) + Bm = 0. 



(8) 



(9) 



3 



The determination of the MR ratio of the system consists of the following three steps. 
First we determine the charge current density je from the continuity equation V-je = 0. Next 
we solve the Eq. for a given local magnetization M in the domain wall and determine the 
spin accumulation m. Finally we determine the voltage drop of the system by integrating 
the electric field E = (jg + 2D • Vm)/(2Co). Since we apply the constant electric current 
density, the voltage drop is proportional to the total resistance, from which we obtain the 
MR ratio of the system. 

Now we move onto the analysis of the MR ratio of the CCP spin-valve with a domain 
wall. The system we consider is schematically shown in Fig. [T](a). A domain wall is confined 
in a nanocontact, which is sandwiched between pinned layer and free layer. We adopt Bloch 
domain wall as a magnetic structure in the nanocontact. The local magnetization M is given 
by M = {cos{7Tz/zh),sm{TTz/zh),0), where z^ is the thickness of the nanocontact (=2nm) 
and z is the height from the bottom of the nanocontact as is shown in Fig. [1] (b). 

The cross-sectional view of the system is shown in Fig. [T] (c). The curve of the contact 
region is modeled by the half ellipse, which is tangent to the ferromagnetic layer. The length 
of the major axis of the ellipse, which is the same as the thickness of the contact, is taken 
to be 2 nm. The length of the semi-minor axis of the ellipse is 1 nm. The pinned and free 
layers are assumed to be 40nmx40nmx4nm rectangles. We also consider the Pt capping 
layer with the thickness 50nm. We set the value of the resistivity of Pt layer as large as 
Iflum so that the total RA without CCP structure is O.lfium^ which is the value reported 

riQ 

in the experiments y| • 

The spin accumulation m is determined by numerically solving Eq. by use of the 
finite element method Q], [sl. The system is divided into hexahedral elements and the total 
number of elements is of the order of 10^. We assume that the free and pinned layers and 
the nanocontact are all made of CoFe. We use material parameters for the conventional 
CPP-GMR spin-valve system; (3 = 0.65, pcoFc = 150f2nm, AcoFc = 12 nm. Apt = 14 nm, 
J = 0.5 eV, where A denotes the spin diffusion length. Especially we adopt the value 
of resistivity of CoFe in the contact region as IQ^m, which is the value reported in the 
experiments [l, 2]. This enhancement of the resistivity in the contact region is due to the 
decrease of the purity of CoFe in the oxidation process. The diffusion constant Dq and the 
relaxation time r^j are determined by the relation Cq = NpDo and A = a/2ts/Do(1 — PP'), 
where we set the density of states at the Fermi level as Np = 7.5nm~^eV~^. 
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In order to examine the MR enhancement effect by the CCP structure, we first vary the 
contact diameter d with the fixed size of the electrode D = 40nm. The dependence of MR 
on d is shown in Fig. [2] (a). The MR enhancement effect by the CCP structure is confirmed 
by the simulation as is shown in the figure, where the MR ratio increases with decreasing the 
contact diameter d from d = 20nm to d = 2nm. If we completely confine the contact region 
(i — 0, no electric current can flow and MR ratio becomes zero. This implies the existence of 
the contact diameter which maximize the MR ratio. Actually we observe sharp peak around 
d ~2 nm in Fig. [2] (a). In practice, in the region where the contact diameter is less than the 
order of 1 nm, the diffusion equation loses its validity and the quantum mechanical approach 
is necessary. In Fig. [2] (b) we replot the MR ratio against the resistance area product RA. 
The value of RA without CCP structure is set at O.lQfim^. It is observed that if we decrease 
the contact diameter, the MR increases with increasing RA and takes its maximum value 
around RA~ IQ^m^. 

Next we vary the size of the electrode D with the fixed contact diameter d=2nm. Since 
enlarging the size of the electrode D corresponds to increasing the distance between nanocon- 
tacts in the NOL, the effect of changing the density of nanocontacts can be realized in our 
simulation by varying the size of the electrode. The dependence of MR on D is shown in Fig. 
[3] (a). The MR is a monotonic increase function of D and is saturated around D ~50nm 
as is shown in the figure. This means that the MR monotonically increases with decreasing 
the density of nanocontacts and is saturated at a certain value of density. In Fig. [3] (b) 
we replot the MR ratio against the resistance area product RA. Since RA increases with 
decreasing the density of nanocontacts, we observe a monotonic increasing dependence of 
MR on RA unlike the case with varying the contact radius d in Fig. [2] (b) where MR has 
maximum value at RA~ IQfim^. Although the MR shows monotonic increasment with RA, 
it is saturated around RA= 1 ~ 2f2/im^. 

In summary, we studied the current-perpendicular-to-plane magnetoresistance of a do- 
main wall confined in a CCP structure made of an NOL. We solved the continuity equations 
for charge and spin currents with the three-dimensional CCP geometry by use of finite ele- 
ment method. We have confirmed that the MR enhancement effect by the CCP structure 
from the calculation of MR with varying the contact radius d and the electrode size D. We 
observed that the MR is maximized at d ~2nm. In our simulation with varying the electrode 
size D, it is observed that the MR increases with decreasing the density of nanocontact and 
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is saturated at a certain value of density. 
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FIG. 1: (a) A CCP spin-valve with a domain wall is schematically shown. A domain wall is confined 
in a nanocontact. (b) The close-up of the contact region is schematically shown. We adopt Bloch 
domain wall as a magnetic structure in the nanocontact M = (cos(7rz/2;^), sin(7r2;/z/i), 0). (c) The 
cross-sectional view of the CCP spin-valve is schematically shown. The curve of the contact region 
is modeled by the half ellipse, which is tangent to the ferromagnetic layers. In our calculation we 
keep the shape of the ellipse fixed and vary the size of the contact by changing the contact diameter 
d. We also vary the size of the electrode D with fixed contact diameter d. 
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FIG. 2: MR ratios of the domain wall CCP spin-valve are plotted with varying the contact 
diameter d. (a) MR is plotted against contact diameter d. (b) MR is plotted against resistance 
area product RA. The value of RA without CCP structure is set at 0.iri//m^. 
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FIG. 3: MR ratios of the domain wall CCP spin-valve are plotted with varying the size of the 
electrode D. (a) MR is plotted against the size of the electrode D. (b) MR is plotted against 
resistance area product RA. The value of RA without CCP structure is set at O.l^lfim^. 
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